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Abstract 

The geometric phase analysis (GPA) algorithm is known as a robust and straightforward technique that can 
be used to measure lattice strains in high resolution transmission electron microscope (TEM) images. It is 
also attractive for analysis of aberration-corrected scanning TEM (ac-STEM) images that resolve every atom 
column, since it uses Eourier transforms and does not require real-space peak detection and assignment to 
appropriate sublattices. Here it is demonstrated that in ac-STEM images of compound materials (i.e. with 
more than one atom per unit cell) an additional phase is present in the Eourier transform. If the structure 
changes from one area to another in the image (e.g. across an interface), the change in this additional phase 
will appear as a strain in conventional GPA, even if there is no lattice strain. Strategies to avoid this pitfall 
are outlined. 
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1. Introduction 

Quantitative measurement of strain on the 
nanometre scale is important for many materials 
studies, and it is increasingly evident that many 
interesting phenomena can only be understood 
through knowledge of the structure at length scales 
on or below the level of individual unit cells. The 
strain associated with an interface is often of spe¬ 
cial interest, both for an understanding of the way 
the materials are produced [1] and for the function¬ 
alities that result [2j [31 01 [5l El [7]. In transmis¬ 
sion electron microscopy (TEM), there are several 
methods that can be used to make such a mea¬ 
surement, including diffraction, holography and the 
analysis of high resolution images that resolve the 
crystal lattice, as compared recently by Beche et 
al. [H]. Eurthermore, in the case of high resolu¬ 
tion TEM images, there are several different ap¬ 
proaches [9], all with origins dating to the 1990s. 
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These can be broadly classified into three differ¬ 
ent types: direct measurement of interatomic dis¬ 
tances in real space mm, extraction of a lat¬ 
tice by comparison to a template mm and anal¬ 
ysis in Eourier space mm- Although they were 
originally applied to conventional high resolution 
TEM images [15] [16], they are increasingly being 
applied to aberration corrected high-angle annular 
dark field scanning mode (HAADE-STEM) images 
[Il[IZl[IHl[l9l[20l[2l]. In principle, the lack of con¬ 
trast reversal with specimen thickness or defocus 
in such images makes them well-suited for strain 
analysis. Eurthermore, atomic resolution measure¬ 
ments of strain and displacement can be combined 
with atomic number (Z) contrast and simultane¬ 
ous spectroscopy such as electron energy loss spec¬ 
troscopy (EELS) or energy dispersive X-ray (EDX) 
analysis to give a wealth of information. Strategies 
to average, or remove, distortions [20] in HAADE- 
STEM images are now available [22] [23] , allowing 
strain analysis to be employed with confidence on 
these data-rich images. Nevertheless, when there is 
more than one atom per unit cell, i.e. in images of 
compound materials, complications can arise that 
can compromise a strain analysis. 
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Here, the implications of atomic resolution on 
the Fourier space technique commonly known as 
geometric phase analysis (GPA) [T3], one on the 
most widely used techniques. GPA is particularly 
attractive for atomic resolution images with sev¬ 
eral atoms per unit cell, since there is no need to 
assign atomic columns and interatomic distances 
to different sublattices, which must be done for 
real space techniques. However, additional precau¬ 
tions must be taken to minimise unwanted effects 
from the microscope and specimen that can pro¬ 
duce incorrect strains, particularly in TEM [IS]. It 
is shown here that additional artefacts often ap¬ 
pear in GPA analyses of atomic resolution images, 
particularly at interfaces of materials of different 
structures, and these can be easily attributed to 
a lattice strain that does not exist. To resolve this 
problem, strategies to measure true strain maps are 
presented. The problem is illustrated using two ex¬ 
amples, an InGaAs-AlAsSb quantum cascade laser 
and an SrRuOs-SrTiOs (SRO-STO) interface, both 
with aberration-corrected ADF-STEM experimen¬ 
tal data and multislice simulations of strain-free 
model structures. 



Figure 1: (a) and (b) show two different example structures 
created as the sum of two peak types, (c) shows the fringes 
from the lowest order Bragg spots of (a) (dotted blue) and 
(b) (solid orange). The highlighted positions of the maxima 
are different and do not correspond to any peaks in the origi¬ 
nal structures, though they are closer to the higher intensity 
peaks. 


2. Theory 

The core principal of GPA mm] is based on the 
concept of a local value Hg (r) of a Fourier compo¬ 
nent associated with lattice fringes described by the 
reciprocal lattice vector g. A digital image with in¬ 
tensities that vary as a function of position, / (r) 
can be expressed by a Fourier series, summed over 
all reciprocal lattice vectors; 

I {r) = ^Hg{r)exTp{2Trig-r), (1) 

9 

where the Fourier components, Hg (r), are given 
by 

Hg (r) = Ag (r) exp (iPg (r)) (2) 

with Ag (r) as the (real) amplitude and Pg (r) 
giving the locally varying phase. The particular 
importance of this phase is that it describes the po¬ 
sition of the lattice fringes and, by comparison to a 
reference lattice fringes, can be used to calculate a 
displacement field and strain in the image. In prac¬ 
tice, the components Hg (fc) are extracted in Fourier 
space by application of a mask (usually Gaussian) 
around the reciprocal lattice point g. An inverse 
Fourier transform yields Hg (r) and the phase of 
the complex image is then calculated as 

Pg (r) = Phase [Hg (r)] - 27rp • r = -27rp • ti, (3) 

where u is the displacement of the lattice fringes 
from the reference planes given by g. Using the 
measured phase, it is then trivial to calculate the 
displacement u in the direction oi g. By repeating 
this process for a non-collinear p, the full displace¬ 
ment field can be measured which is then differen¬ 
tiated to get the strain m- 

An implicit assumption in this process is that the 
change in phase for the lattice fringes g is simply 
due to the displacement of the atoms in the mate¬ 
rial. This is almost always correct when the reso¬ 
lution of the image is only sufficient to resolve the 
crystal lattice. Nevertheless, a high resolution im¬ 
age should be considered a convolution of a lattice 
(a mathematical array of points described by a 2D 
Dirac comb function IH = 6 {r — Pj)^ where Pj are 
vectors describing the lattice points) and a basis 
image / (e.g. the image of a single atom, placed at 
each lattice point): 
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/ (r) = HI (g) /. 


( 4 ) 












Prior to aberration correction, the resolution of 
many TEM or STEM images was sufficiently low 
that the basis image was closely approximated by 
a simple sine wave. However, at higher resolutions 
the basis image / also contributes to P^, i.e. the 
phase obtained by GPA is no longer simply related 
to the lattice displacement, but the convolution of 
lattice and basis. Erom this point of view, it is 
easy to see that if the basis image is represented 
by a locally varying function, /(r), that changes 
from one place to another (e.g. across an interface) 
there will in general be a change in phase - even in 
an image with an unchanging lattice. 

This effect can be understood by considering a 
high resolution image / (r), comprised of a lattice 
convoluted with a basis image that contains two 
atoms of different types. An equivalent description 
is the sum of two images, each consisting of a lattice 
convoluted with a basis image containing only one 
atom. The image can then be described as the sum 
of two sublattices, A and B, each represented by a 
Eourier series: 



Figure 2: (a) ADF-STEM image along the [110] axis of 

Ino.rGao.sAs-AlAsSb layers. The overlay on the right high¬ 
lights the structure and lattice planes in the image, (b) In¬ 
tensity profile from the overlaid regions in (a), following the 
direction of the arrow. 


/ (r) = y] exp {2wig ■ r)+ 

® (5) 

exp { 2 'Kig • r) exp { 2 'Kig • -y), 

9 

where the A lattice has been taken as the origin 
and V is di vector describing the displacement be¬ 
tween the sublattices. Assuming that the images of 
the A atoms are identical to the images of the B 
atoms, only differing in intensity, it can be written 
that 


= aHg ( 6 a) 

Hf = pHg ( 6 b) 

where Hg (r) are the Eourier coefficients of a nor¬ 
malised image of a monatomic crystal. Using this, 
Eq. can then be rewritten as a single Eourier 
series 


I{r) 

where 


Cg exp { 2 'Kig • r) exp {i(j)) 

9 


(7) 


Cg = Hg^Ja^ ^ 2a^ cos {2Kg • v) 
(3 sin {2Kg ■ v) 

ncp = ---T--. 

a-\- p cos {2Kg • v) 


( 8 a) 

( 8 b) 


Note that Eq. [^has exactly the same form as Eq. 
[I] but with an additional factor that is related 
to the basis image. Thus, when applying GPA to a 
high resolution image with multiple sublattices, the 
calculated phase will be 


P 9 (^) = • u - 


(9) 


where there is the term (p (r) describing contri¬ 
butions from the basis image in addition to the lat¬ 
tice strain. Note that if the basis image does not 
change as a function of position, 0 simply adds a 
constant phase across the whole image and has no 
effect on the strain components obtained by differ¬ 


entiating. It is readily apparent from Eq. [ 8 b| that 
there are many possibilities for additional phase 
shifts that can produce unwanted artefacts in strain 
maps. The situation becomes more complicated 
when applying GPA to atomic resolution images 
with several sublattices A, B, C... when all parame¬ 
ters change as a function of position, i.e. intensities 
a{r)^ /3 (r), 7 (r). .. and relative sublattice displace¬ 
ments Vb (^), Vc {r).... Nevertheless, the phase 
shift does not affect all g in the same way and it 
can be possible to choose appropriate g such that a 
true strain map can be obtained. 

Eor a GPA map to show only strain, the addi¬ 
tional phase (j) must be zero. Eor a biatomic unit 
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Figure 3: (a-c) Syy strain measured using the Bragg spots highlighted in the inset FFTs. (d) strain profiles from the boxed 
regions in (a) (dashed black), (b) (solid orange) and (c) (dotted blue). 


cell, inspection of Eq. [^shows that this occurs when 

gv = n, (10) 

where n is an integer (including zero). If this 
is not the case, an additional phase shift will be 
present and can appear as a false ‘strain’ parallel 
to V with a magnitude proportional to g ■ v. In an 
image with two sublattices, described by Eq. it 
is easy to see that when the A sublattice is much 
brighter (a ^ /3), 0 ^ 0 and the phase represents 
the A lattice; conversely when a (3^ (p ^ 2Tvg • v 
and the phase is that of the B lattice. Simply put, 
the sublattice with the larger amplitude has most 
influence on the phase, as demonstrated in Eig. 

3. Experimental 

All experimental images were taken using a dou¬ 
ble CEOS aberration-corrected, Schottky emission 
JEOL ARM-200E microscope operating at 200 kV. 
Aberration corrections were applied to third or¬ 
der and measured to fifth order [251 |26l [27]. In 
STEM mode, a convergence semiangle of 22 mrad 


was used with a JEOL ADE detector with inner and 
outer collection semiangle of 45 and 180 mrad re¬ 
spectively. To eliminate scan distortions, the spec¬ 
imen was left until drift was less than 20 pms“^ 
and a series of up to 50 images was collected, 
each with a short pixel dwell time (< 2ms). A 
high quality image was obtained by aligning sub¬ 
sequent images using normalised cross correlations 
and then summing the complete series. GPA was 
performed using a program developed in-house; all 
basis directions for the strains are set the same 
as the image bases. Simulations were performed 
using clTEM, an open source, GPU accelerated 
multislice program m- The simulation parame¬ 
ters were matched to the microscope aberrations 
as measured by the GEOS control software; ther¬ 
mal diffuse scattering was modelled using the frozen 
phonon method with 15 configurations m- 
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Figure 4: (a) Simulation of Ino.rGao.sAs-AlAsSb interfaces, (b-d) £yy strain measured using the Bragg spots highlighted in 
the inset FFTs. 

4. Results and discussion 

4.1. InGaAs-AlAsSb interfaces 

An ADF-STEM image of two interfaces in an epi¬ 
taxial layer structure for a quantum cascade laser 
is shown in Fig. [^a). In this image, taken along 
the [110] crystal axis, the group III and group V 
atom columns are readily resolved as a ‘dumbbell’, 
with the lower part corresponding to the group III 
sublattice (A) and the upper part corresponding to 
the group V sublattice (B). The displacement of 
the two sublattices in this projection is t; = | [ 001 ] 

(there is also a displacement along the incident elec¬ 
tron beam direction, which has no effect on the im¬ 
age). The dark, horizontal, central band is a thin 
AlAso. 8 Sbo .2 layer, roughly 3 monolayers in thick¬ 
ness, between Ino. 7 Gao. 3 As layers (top and bot¬ 
tom). In such images one is often interested in the 
ability to achieve the epitaxial design, and strain 
measurement can be a crucial part of this assess¬ 
ment. In this image, the relative intensities of the 
sublattice A and B swap as shown in Fig. [^b). The 
occupancy of the A sublattice switches from In-fGa 


(brighter) to A1 (fainter), while the B sublattice 
changes from As (fainter) to As+Sb (brighter). 

Figure has the characteristics of an image that 
will contain phase shifts in Fourier components that 
are caused by the basis image, rather than the lat¬ 
tice. Figure [^a) shows the Syy output of GPA per¬ 
formed using Fig. |^a), using the g = 002 and 
g = 022 spots. Since g v = ^ for g = 002 , this com¬ 
ponent is affected by the changes in the basis image. 
There is an apparent strain of 10% at the interfaces 
that does not agree with a visual inspection of the 
image. Gonversely, g • v = 0 for the g = 220 and 
it is insensitive to this effect. Naively, one might 
hope that this issue may be overcome by avoiding 
the g = 002 component and using, for example, the 
g = 111-type Fourier components. Nevertheless, 
in this case g v = Therefore the g = 111 and 
= 111 spots are also expected to produce specious 
strains at the interfaces, as is the case in Fig. |^b). 
The effect is much smaller since g ■ v is small, and 
may be easily overlooked. A true strain map (Fig. 
I^c)) is only obtained by choosing p-vectors that 
obey Eq. 10 for example g = 004 and g = 220. 
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Figure 5: (a) ADF-STEM image of an STO-SRO interface, 
as highlighted by the structure on the right along with the 
lattice planes, (b) Intensity profile of the boxed regions in 
(a) following the overlaid arrow. 


This shows the true strain at the interface to be 
only r\j 1 %. 

It is apparent in Fig. that there is some 
real lattice strain that is combined with the arte¬ 
fact produced by changes in basis image. In or¬ 
der to demonstrate phase shifts without any lat¬ 
tice strain, a similar procedure was performed on a 
multislice-simulated ADF-STEM image of a strain- 
free Ino. 7 Gao. 3 As AlAso. 8 Sbo .2 model heterostruc¬ 
ture (Fig. 1^. The GPA maps use the same sets of 
p-vectors as used in Fig. The apparent strains at 
the interfaces can be seen to have a similar sign and 
magnitude as in the experimental images of Fig. 
although they are more visible (mainly because the 
interface is perfectly abrupt) and are zero in Fig. 

m 

4.2. STO-SRO interface 

ABO 3 perovskites are the focus of intense study 
using aberration-corrected STEM [siEniisi], and as 
a compound materials they also exhibit erroneous 
strain when examined using GPA. This has on occa¬ 
sion been interpreted as a real interfacial strain [18] . 
Most investigations to date have used ADE-STEM, 
which is insensitive to the oxygen atoms, and so can 
be considered to produce biatomic images. How¬ 
ever, annular bright field (ABE-STEM) is increas¬ 
ingly popular [32j [33l |3l| since these images also 


show oxygen atoms [35ll36]. In (100) ABE-STEM 
images of such materials there are effectively four 
sublattices, meaning that great care must be taken 
in strain analysis. The analysis here is restricted to 
the simple case of a [001] ADE-STEM image, with 
only two sublattices related hy v = ^ [ 110 ] (note 
this choice of v is somewhat arbitrary since it can 
be defined modulo any lattice vector). An interface 
between STO (bottom) and SRO (top) is shown in 
Fig. i a) and, although the Sr sublattice remains 
constant across the interface, the A and B sublat¬ 
tices show a form of inversion in the basis image as 
Zxi < Zsr in STO while in SRO Zru > Zsr- 

In this case it is possible to produce artefacts that 
appear as large shear strains (that are clearly un¬ 
physical) as well as those that appear perpendicu¬ 
lar to the interface. Eigurej^a) shows the axial and 
shear GPA maps produced using the g = 100 and 
g = 010 Eourier components, both of which give 
g • V = Since this gives a basis image-induced 
phase shift of tt, strong artefacts in a GPA strain 
map are to be expected. Indeed, at the interfaces in 
Eig. [^a), strains greater than 20% are found both 
in the Syy and Syx components. As before, such 
artefacts can be avoided by choosing spots such as 
g = 110 and g = 110 that give g ■ v = 1 and zero 
respectively [^b). It is evident that there is in fact 
minimal strain in the image (< 2 %). Analysis of 
a simulated image from a strain-free model struc¬ 
ture produces the same results (Eig. with strains 
measured as zero in Eig. [^c). 

The strain artefacts shown in Eigs. 4-7 have 
been chosen to be large and easy to distinguish; 
the changes in basis image are large and abrupt, 
while g vectors were generally chosen to maximise 
the additional phase shift (j). However this is not 
always the case and the effect may be quite subtle. 
In particular, it is important to note that the effect 
will still be present even in cases where the relative 
intensities of the sublattices do not reverse, Eor a 
biatomic unit cell, it is quite clear that any change 
in (X or /3 in Eq. will change 0, and this will 
appear in a GPA map as a strain if inappropriate g 
vectors are chosen. 

It should be noted that the artefacts in strain 
maps are as a result of extra information that need 
not always be disregarded. Eor example, the rate 
at which the basis image changes will provide a 
phase shift that can be used to indicate the sharp¬ 
ness of an interface. Alternatively, an indication of 
sub unit-cell displacements (e.g. in ferroelectrics) 
across a domain wall can be obtained, even if there 
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Figure 6: (a), (b) GPA strain maps for the Sxx components and Syx components (lower insets) produced using the Bragg spots 
highlighted in the upper inset. 



Figure 7: (a) Simulation of an STO-SRO interface, (b), (c) GPA strain maps from (a), showing the Sxx components and Syx 
components (lower insets). The inset FFTs show the Bragg spots used for the analysis. 


is no lattice strain. With careful consideration of 
what the Fourier components represent, it is pos¬ 
sible to exploit this extra information, rather than 
be inhibited by it. 

5. Summary 

It has been shown that in atomically resolved im¬ 
ages, one must take care to distinguish changes in 
the lattice from those in the basis image. Lattice 
strain, changes in relative intensity of sublattices 
and sub unit cell displacements all produce phase 
shifts in Fourier components. In several cases rep¬ 
resented here the phase shifts are large and should 
be obvious as artefacts (e.g. a 20% shear strain con¬ 
fined to sub-nm dimensions should be perceivable 
by eye). However, in others they are quite sub¬ 
tle and could easily be mistaken for real interfacial 
strains (e.g. analysis of (110) III-V semiconductor 
images using 111-type g vectors). In images with 
two sublattices, the problem can be avoided by only 


using p-vectors which satisfy the rule g where n 
is an integer. 
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